Magnetoelastic-magnetoelectric phase transitions in multiferroic BiFeC>3 
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Our measured dielectric constant and mechanical response of multiferroic BiFeOa indicate four 
phase transitions below room temperature. Features correlate with those reported at 50K (from a 
peak in the zero-field-cooled magnetic susceptibility) and 230K (from splitting between field-cooled 
and zero- field-cooled magnetic data [T]), and those at 140K and 200K (from magnon light scatter- 
ing cross sections [5|)- The primary order parameter is not the polarization in any of the low-T 
transitions. Instead, the transition near 200 K shows strong elastic coupling, while that at 50K is 
fundamentally magnetic, but magnetostrictively coupled to the the lattice. The low-T phase tran- 
sitions display glassy behaviour. A further anomaly at 140K interpreted as spin reorientation [2] [3] 
shows only weakly in dielectric and mechanical studies, indicating that it is predominantly magnetic 
with little coupling to any of the other order parameters. 



PACS numbers: 75.80.+q; 77.90.+k; 77.80. Bh;81.40.Jj 

Magnetoelectric multiferroics are materials where fer- 
roelectricity and magnetism coexist and are coupled. 
These materials are attracting increasing interest due 
to their possible applications as magnetic detectors and 
multi-state memory devices, as well as for the rich physics 
behind the fundamental mechanisms for magnetoelectric 
coupling. Among the materials being most actively stud- 
ied, BiFeOs (BFO) stands out for being the only known 
room T magnetoelectric multiferroic (ferroelectric, an- 
tiferromagnetic) , and has some very attractive features 
such as huge remnant polarization (around 100 ^C/cm 2 ) 
in both single-crystal |U [5] and thin- film form [5) [7] and 
coupling between the directions of the ferroelectric po- 
larization and the sublattice magnetization |5]. For this 
reason, BFO has been touted as a possible active mate- 
rial in the next generation of lead-free ferroelectric ma- 
terials [5]] and also as the active layer in magnetoelectric 
memories based on exchange bias [10] . 

BiFeOs has a rhombohedrically distorted perovskite 
structure (space group R3c) at room T [11] which is 
both ferroelectric (with polarization along the (111) di- 
rections) and antiferromagnetic. The antiferromagnetic 
symmetry is itself very exotic, with a long period (62 nm) 
cycloidal helix of spins, inconmensurate with the lattice 
periodicity [12j [13]. The cycloidal rotation cancels al- 
most all macroscopic magnetization, and hence BFO is 
antiferromagnetic; upon heating, it becomes paramag- 
netic at Tat=640K [H]. Around HOOK BFO shows a 
transition from rhombohedral to orthorhombic symme- 
try [7] or perhaps monoclinic [8], and around 1200K it 
becomes cubic. It is not yet certain which of the two 
phase transitions is the ferroelectric-paraelectric one. On 
the other hand, there is evidence that the transition to 
the cubic phase is accompanied by the onset of a metal- 
insulator transition [T3 US]- While the high T phase 
diagram appears to be quite rich, there has been compar- 
atively little reported for the low-T properties, although 
it has been suggested [T7| [T5] that spin-glass behaviour 



occurs amorphous BFO below 200K. However, three re- 
cent works [TJ HJ [3] indicate further complexity in the 
low-T nature of BFO. 

The first of these [1] reports the magnetic proper- 
ties of BFO single crystals at low temperatures. It was 
found that the field-cooled and zero-field cooled magne- 
tization became different at temperatures below w 250K, 
as can subsequently be seen in earlier data [HE |2"U] , 
Furthermore, there is a distinct peak in the zero-ficld- 
cooled magnetization around 50K. Below this lower T, 
the temperature-frequency dependence of magnetic sus- 
ceptibility indicates a possible spin glass transition, al- 
though the authors noticed that the behavior was un- 
usual in that susceptibility increases with increasing fre- 
quency. It was also reported that the critical exponent 
describing the slowing down of the glassy dynamics {zv), 
was very small (« 1 .4) , much closer to the value expected 
in a mean-field system (where zv=2) than in a classic 
Ising-type magnetic spin glass (zv—l-lQ) [2T]. On this 
basis, it was suggested that the spin glass transition may 
be coupled to a long-range order parameter responsible 
for its mean-field behaviour. 

Later, two contemporaneous and independent Ra- 
man studies reported electromagnon Raman spectra in 
BFO [2, 3 . Singh ct al. |2J noted two anomalies in the 
Raman frequency, intensity, and linewidth: one strong 
at «200K, one very weak at «140K. Cazayous et al. re- 
ported a strong anomaly at 140K, and only a keen eye 
might discern a very small anomaly at 200K. It is in- 
teresting that the strength of the anomalies was so dif- 
ferent. However, in more recent work they found both 
transitions with equal magnon cross section divergences 
by changing their scattering geometry, showing that the 
differences are not sample-dependent but orientation- 
dependent. Two separate spin reorientation transition 
temperatures might be expected since the spins begin to 
rotate out of a plane at the upper one, while at the lower 
one the spins have rotated fully 90°, so that they are or- 
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FIG. 1: Frequency dependence of the dynamic mechanical 
modulus and tan<5 of BiFe03 as a function of T. The large 
peak in tan<5 and concomitant relaxation in modulus on heat- 
ing through 200K shows significant frequency-temperature 
dispersion. 



thogonal to the original plane. Singh et al.'s [TJ data are 
compatible (in respect of both peak intensity divergence 
and linewidth narrowing) with the original study of spin 
fluctuations in uniaxial antiferromagnets by Schulhof et 
al. [22]. They found that fluctuations along the uniaxial 
direction diverge and exhibit critical slowing down (spec- 
tral narrowing) approaching T c , with critical exponents v 
= 0.63 and 7 (susceptibility) — 1.24 for the longitudinal 
fluctuations and v — 0.63 and 7 = 1.47 for the trans- 
verse fluctuations. Note that the transverse fluctuations 
are not truly divergent but extrapolate to a divergence 
w2K below T c . 

In order to understand further the nature of these low- 
T anomalies, we have studied them from the point of 
view of the two other order parameters of this multi- 
ferroic material: ferroelasticity and ferroelectricity. Dy- 
namical Mechanical Analysis of ceramic BFO bars was 
performed in a three-point bending configuration using 
samples that were 2mm wide, 0.5mm thick, and 5mm 
between knife edges. Force between 30 and 630 mN were 
applied dynamically, generating strains of the order of 
10~ 4 , at frequencies between 0.6 and 18 Hz and Ts be- 
tween 120 and 290K (Fig. [IJ. The Youngs modulus and 
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FIG. 2: The mechanical modulus and tan<5 of BiFe03 mea- 
sured at 2Hz. A small anomaly in tan<5 can be seen at 140K. 



loss angle (indicative of a mechanical dissipative process, 
such as domain wall friction [23J show stiffening of «10 
% upon cooling accompanied by a simultaneous peak in 
dissipation, tan5, indicative of transition from a dynam- 
ically relaxed system above 230K to an unrclaxed state 
below this temperature. There is strong frequency de- 
pendence of the mechanical relaxation T, signaling ei- 
ther a true glassy transition or else a defect-related re- 
laxation. The limited range of frequencies accessible in 
our measurement makes it difficult to compare Vogel- 
Fulcher (glassy freezing) and Arrhenius (thermal activa- 
tion of defect states) models for the response. None the 
less, analysis of the frequency dependence of the peak in 
tand around 200K is consistent with an activation energy 
of 0.59(9) eV which is smaller than that expected for me- 
chanical relaxation due to oxygen vacancies in perovskite 
oxides (typically around 0.7 to 1.1 eV [2"4"]1. 

The DMA data also show a very weak anomaly, seen 
most clearly in the data for tan5 collected at 2Hz (Fig. |2| , 
at 140K. The weakness of the 140K feature and strength 
of the one around 200K contrasts markedly with the rel- 
ative scale of the features seen at these Ts in the elec- 
tromagnon spectra, which may indicate variable coupling 
between strain and the primary order parameter at the 
lower-T anomaly. 

The dielectric constant was also measured as a func- 
tion of T and frequency for single crystals and for ceram- 
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FIG. 3: Dielectric constant of BiFeOs at low-T. A small 
anomaly at 140K can be seen in the single crystal data at the 
highest frequency. 

ics (see Fig. [3]) . The single crystals were cut with faces 
perpendicular to {110} pseudocubic planes. The dielec- 
tric constant shows subtle but unambiguous anomalies 
around 50K, 200 and 230K. The fact that they are subtle 
may explain why they may have gone unnoticed in spite 
of extensive research on this compound in recent years. It 
also indicates that the primary order parameter in these 
transitions is almost certainly not the polarization, and 
that they reflect indirect coupling to the primary order 
parameter. While a single broad peak, strongly depen- 
dent on frequency, is seen around 200K in DMA, two 
resolved anomalies (independent of frequency) are seen 
in the dielectirc data. 

We note that there is no strong evidence of the 
anomaly at 140K except at the highest frequency used 
(400 kHz). This suggests that further studies in the 
MHz regime should be carried out on single crystal sam- 
ples; it may be that this transition is insensitive to low 
frequency probes. The dielectric constant (capacitance) 
shows some low-frequency dispersion, probably due to 
space-charge, but nevertheless the position of the anoma- 
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FIG. 4: Dielectric constant and loss of BiFe03 at low-T 1 , 
measured at 10kHz. 

lies appears to be relatively frequency-independent, con- 
trary to the glassiness observed in the magnetic and me- 
chanical data. It is also worth mentioning that although 
the space charge increases the dispersion and the diel- 
electric losses, these remain low (tan<5<0.1) (Fig. [4} and 
furthermore (perhaps surprisingly) there is no obvious 
sign of change in the losses at the critical temperatures, 
ruling out a conductive artifact |25] . 

Both the mechanical and dielectric data show only 
very weak evidence of the transition at 140K, which sug- 
gests that this transition is predominantly magnetic. The 
transitions above «200K, on the other hand, have now 
been observed in different samples (single crystals and 
ceramics) and with five different techniques (magnetiza- 
tion, susceptibility, Raman, impedance and mechanical 
spectroscopy). If either is due to defects, they are of 
a very pervasive but undetectable nature: XRD studies 
did not detect any impurities or second-phases in any 
of our samples. Furthermore, the mechanical response 
seems very large to be caused by to defects alone, sug- 
gesting instead that the primary order parameter may 
be a structural one, coupled weakly to the polarization 
(hence the dielectric anomaly) and to the magnetization 
(hence the electromagnon and magnetization anomalies). 
In the dielectric data two small discrete anomalies at 200 
and 230K are seen while the mechanical relaxation at this 
temperature is broad and large. 

The low-T anomaly indielectric response at 50K is 
weak but nonetheless clear, showing that the magnetic 
order parameter is coupled to the polar one. Since the 
latter is a long-range order parameter, this coupling may 
account for the apparent (and very unusual) mean-field 
nature of the low-T spin glass transition. We note also 
that if this material is indeed a spin-glass, it would be 
acentric. Fischer and Hertz (5^1 have emphasized that 
no published theories apply to spin glasses lacking an 
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inversion center, and further, that such glasses cannot 
possibly be Ising-like. 

Our results indicate that BiFeC>3 is remarkably sim- 
ilar to the orthoferrite ErFeC>3, which has T/v = 633K 
and spin-reorientation transitions at 90K and 103K [com- 
pared with T/v = 643K and spin reorientation transitions 
at WOK and 200K in BiFe0 3 . According to the the- 
ory of Koshizuka and Ushioda [27] one might expect the 
magnon frequency / to drop to zero at T reorientation and 
the magnon integrated cross section to diverge due to en- 
hanced thermal population kT/h/; however Koshizuka 
and Ushioda report [PRB 22, 5394 (1980)] a smaller 
(50%) magnon frequency decrease and consequently a 
cross section increase of only about x5, and attribute this 
to phonon-magnon coupling (magnetostriction). Singh 
et al.'s 2 magnon frequency decrease is only «5 % and 
the peak scattering amplitude increases by a maximum 
of x5 to xlO - but the linewidth decreases by at least 
x5 (resolution limited). Thus the integrated cross sec- 
tion changes by only a small amount. This is compatible 
with the strong mechanical losses we see around 200K 
(and somewhat weaker at 140K) and implies very large 
striction. This is compatible with respect to linewidth 
narrowing and peak amplitude divergence with the spin- 
fluctuation models of Koshizuka and Ushioda and Schul- 
hof et al. [251 [57]; however, the linewidth narrowing and 
peak amplitude divergence imply critical fluctuations not 
explicit in their original papers. 

In summary, we show that the anomalies in bismuth 
ferrite at cryogenic Ts can be interpreted as phase tran- 
sitions at 50K (magnetic, but glassy and with magne- 
toelectric coupling), WOK (dominantly magnetic), 200K 
(magnetoelastic), and 230K (magnetic; glassy). The cou- 
pling between magnetic, electric and elastic order param- 
eters varies according to the length scale of the transition 
dynamics in each case: dipolar and strain coupling leads 
to the possibility of long-range correlations dominating 
the spin glass behavior, resulting in a rare example of an 
acentric mean-field spin glass in BFO. It would appear 
that a new model for non-Ising spin glasses must now 
be developed if the low-T magnetoelectric properties of 
BiFeC>3 are to be adequately explained. 

We are grateful to R Palai and H Schmid who provided 
ceramic and single crystal samples of BFO. 
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